Wetting on Random Roughness: the Ubiquity of Wenzel Prewetting 
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The wetting properties of solid substrates with macroscopic random roughness are considered as 
a function of the microscopic contact angle of the wetting liquid and its partial pressure in the 
surrounding gas phase. It is shown that Wenzel prewetting, which has been recently predicted for 
a rather wide class of roughness profiles derived from Gaussian random processes by a general dis- 
tortion procedure [30] , should in fact be ubiquitous and prevail under even much milder conditions. 
The well-known transition occurring at Wenzel's angle is accompanied by a prewetting transition, 
at which a jump in the adsorbed liquid volume occurs. This should be present on most surfaces 
bearing homogeneous, isotropic random roughness. 

PACS numbers: 68.05.-n; 68.08.-p; 05.40.-a; 64.75.-g 



While the physics of wetting and spreading on ideally 
smooth surfaces has meanwhile reached a status of ma- 
ture textbook knowledge [JUI], the wetting properties of 
randomly rough solid substrates are still poorly under- 
stood. This is in part due to the vast range of scales to 
be covered, which extends from the nanometer scale of 
atomic roughness to the scale of millimeters, just before 
gravity comes into play. There has already been a lot of 
work concentrating on small scales, where the interplay 
between interface topography and van der Waals forces 
play a role [5HTT] , but the larger scale, ranging from seve- 
lar microns to about a millimeter, which corresponds to 
typical roughness encountered in practical situations, has 
been dealt with so far only for a small class of rather spe- 
cial cases. Most authors have tried to model rough sur- 
faces as Gaussian random processes [T2HPT)] , or to capture 
single aspects of wetting using simplified model geome- 
tries [DM]. 

Since the most common reasons for roughness, like 
wear and weathering, or even deliberate ones like etch- 
ing or sand blasting, are the effect of very many more 
or less independent local attacks, the application of the 
Gaussian random model appears at first glance as the 
most natural, and has consequently been used in a large 
number of studies [T2rtT51 1251 125] . However, it has been 
recently shown that even slight deviations from purely 
Gaussian roughness give rise to qualitative changes in the 
wetting behaviour of the substrate [30] . Using a simple 
model of a distorted Gaussian distribution, it was shown 
that a wetting phase transition appears for a large class of 
systems which is absent for undistorted Gaussian rough- 
ness. This transition, which was termed Wenzel prewet- 
ting, may be of large potential interest in many areas of 
technology, as well as geosciences and biology. The gener- 
ality of that study, however, was quite limited, since the 
underlying random process was still assumed Gaussian, 
and ample use of multivariate normality had to be made 
[T21 [TBI 1251 [25] , In the present paper, it is shown that 
Wenzel prewetting is in fact a universal feature of ran- 
domly rough surfaces, under only very mild conditions 
which arc fulfilled by most customary surfaces. 

Wenzel was the first to report on systematic studies 



of wetting on randomly rough surfaces [31] . He charac- 
terized the roughness by a single parameter, r, which he 
defined as the ratio of the total substrate area divided 
by the projected area. Obviously, r > 1, and r = 1 cor- 
responds to a perfectly smooth surface. The free energy 
which is gained per unit area when the rough substrate is 
covered with a liquid is then given by r(7 Sff — Jsl)> where 
jsi and "f sg are the solid-liquid and solid-gas interfacial 
tension, respectively. If this gain is larger than the sur- 
face tension of the liquid, 7, we expect a vanishing macro- 
scopic contact angle, because covering the substrate with 
the liquid releases more energy than is required for the 
formation of a free liquid surface of the same (projected) 
area. More specifically, force balance at the three-phase 
contact line yields 
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where 9 = arccos [(7 sg — 7 s z)/7] is the microscopic con- 
tact angle [32] • When 9 is reduced to 9w = arccos(l/r), 
9 ma cro vanishes, and the substrate is covered with an 
'infinitely' thick liquid film. In the present paper, we 
will discuss this transition in some detail, both at liquid- 
vapour coexistence and below the saturated vapour pres- 
sure. 

We describe the rough solid substrate by a function 
/(x,y), which shall approximate the actual physical sur- 
face at any required precision, but be mathematically 
smooth, such that V/ and A/ exist everywhere. The 
roughness is assumed to be homogeneous and isotropic, 
i.e., its statistical parameters shall be the same every- 
where on the sample, and independent of rotation of the 
sample about its normal axis. A small amount of liq- 
uid deposited on this substrate will make an interface 
with the surrounding gas, which is described by a second 
function, g(x,y). The support of g is the wetted area, 
which we call W. Continuity of the liquid surface as- 
sures g = f on the boundary of W, i.e., the projection of 
the three-phase contact line, henceforth denoted by dW 
(cf. Fig.gi). 

As the amplitude of most natural roughness is much 
smaller than its dominant lateral length scale, we assume 





FIG. 1: (a) Top view of the sample, showing the wetted areas 
in grey, the bare substrate in white. The normal vector to 
<9W, n, lies in the (x, j/)-plane. (b) Two sets of contour lines 
of / at heights h and h + dh. The hatched area between the 
lines is equal to p(h) dh. 



that 



|V/| « 1 



(2) 



which allows for substantial simplifications. The same 
shall hold for g. The contact angle with the substrate, 9, 
yields the boundary condition 



|V( ff -/)| wtant 



(3) 



which is to be fulfilled everywhere on <9W, to first order 
in 9, V/, and Vg. Applying Green's theorem to (g — /), 
we obtain 



n • V(. 9 - /) ds 



dW 



A(g - /) d 2 x (4) 



w 



where s is the distance along <9W, n its unit normal vec- 
tor, and x = (x,y). Since g = f on dW, V(a — /) is 
everywhere perpendicular to dW. Hence eq. pi) may be 
written as n • V(g — /) ss 6, and eq. Q can be recast into 



[2H - A/] d 2 x = 



(5) 



w 



in which I denotes the length of dW, and H sa \Ag is 
the mean curvature of the liquid-vapour interface. The 
latter is given by the Kelvin equation, 



H = 



ksT p s 
In — 

2<yv m p 



(6) 



where p is the partial pressure of the adsorbed liquid 
in the surrounding gas phase, p s is its saturated vapor 
pressure, v m its molecular volume, and ks is Boltzmann's 
constant. It is thus a convenient measure for the devia- 
tion from liquid-vapour coexistence. 

Before we can exploit eq. ([5]), we derive a few useful 
relations. Let p(f) is the height distribution of /(x), with 
normalization J p(h)dh — \S\ (total sample area). Then 
the total area between the contour lines at / — h and at 



/ = h + dh is given by p(h) dh, which corresponds to the 
hatched area in Fig. [T|b. The average slope on that set, 
ai (h) = (\Vf\) h , is given by 



cri 



/|V/W|dr 
p(h)dh 



(7) 



where dr = dsd/i/|V/(/i)| is the differential of the 
hatched area. Hence for the total length of the contour 
line at height h we obtain 



L(h) = / ds = ai(h)p(h) 



(8) 



In order to express the integral over A/ which appears 
in eq. ([5]), we apply Green's theorem again, this time to 
the area enclosed by a contour line. This yields 



C(h) 



A/d 2 x= / n-V/ds 



DC 



(9) 



where C(h) denotes the set {x | /(x) < h}, and dC 
its boundary, i.e., the contour line itself. Introducing 
C2(h) = (|V/| 2 )/j, we readily see that eq. (J9J) can be 
rewritten as 



C(h) 



A/ d 2 x = a 2 (h) P {h) 



(10) 



In order to fulfill the boundary condition, eq. |3]), the 
vertical position of the three-phase contact line, which 
may be symbolically written as f(dW), will vary along 
dW about an average value, (f(dW)). The projection of 
the contact line onto the plane will thus approximately 
follow the contour line at /(x) = (/(<9W)), with ex- 
cursions towards both the outside and the inside of W. 
These will in cases represent detours, sometimes short- 
cuts with respect to dW. As a reasonable approximation, 
we may thus use I ps L((f(dW))) for the length of the 
three-phase contact line. Similarly, we set 



d 2 . 



)V 



P (f) d/ = W(h) 



(11) 



for the wetted sample area, with h — (f{dW)). Inserting 
these expressions in eq. (p™, we obtain 



2HW(h) 



<Ji{h) 



L(h) 



(12) 



This allows, if p(h), <7i(/i), and 02 (ft) are known from ex- 
perimental characterization of the sample, to determine 
the adsorbed amount of liquid as a function of 9 and H. 
Up to this point, we have not made any specific as- 
sumption about the roughness profile, except its being 
sufficiently shallow for the approximations made above 
to hold. Now we shall go one step further, observing 
that roughness profiles generated by wear, weathering, 
erosion, etching, sand blasting, or similar processes will 



invariably have a finite codomain (we disregard fissures 
and cracks here). In other words, the support of p(h) is 
the interval [h_, h + ] Vx, where h_ represents the depth 
of the deepest trough, and h+ the height of the highest 
elevation on the sample. This has severe consequences 
of o"i and 02, as both must go to zero as h — > h±. 
To see that, consider the distribution of minima, p,_(h). 
Clearly, /i_ — > for h —t /i_ , and we may choose to write 
p~(h) ~ l^o{h — hJ) v ', with some v > 0. Since each min- 
imum is parabolic to first order, one can readily verify 
that it contributes Ait(h — h ) to o\, if h is the depth of 
the minimum. We thus have 



o\(h) < 47r / p_(z)(h- z) dz 



Mo(/i-fe-r+ 2 
+ l)(i/ + 2) 



(13) 



The '<' sign has been used because there may be saddle 
points (or even maxima) occurring at elevations between 
/i_ and h, which can only reduce o\. As eq. ( 13) shows, 



<7i (h) is bounded from above by a function which vanishes 
at least quadratically as h — > h- . An analogous result is 
obtained for h — > h+. 

Next we consider the ratio 01/02- With the abbrevia- 
tion p = |V/|, we have 



01(h) = 2n p 2 q(h,p)dp 



(14) 



where q(h, p) is the distribution of slopes, sampled at 
height h. Similarly, we have 



o 2 (h) = 2ir / p 3 q(h,p) dp 



As a consequence, 

12(h) 
<ri(h) ' 



J p 3 q(h, p) dp 
J p 2 q(h,p)dp 



"^ Pmax I," J 



(15) 



(16) 



where p m ax is the maximum slope encountered at eleva- 
tion h. Since this vanishes as h — > h- as o\ does, so will 



the ratio 02/01, as eq. (16) shows 



Aside from these global properties, both o\ and 02 are 
expected to be largely featureless, due to the general fact 
that the processes leading to roughness exhibit only very 
limited lateral correlation. For any pronounced feature to 
develop in Oi , distant places on the sample would have to 
'conspire' to contribute to that feature at the same depth. 
This can happen only for composite surfaces, where the 
roughness topography may penetrate through a coating 
or other stratigraphic variation of material properties. 
Such ramifications are interesting to consider for practi- 
cal purposes, but well beyond the scope of the present 
article. 

The generic shape of the function 



A(h) = 



02(h) 
01(h) 



(17) 



which appears in eq. (12) is sketched in F ig. |2| accord- 
ing to the discussion above. Following eq. (|12[), the film 
thickness at coexistence (H — 0) can be derived from 
the zeros of A, of which there are either two or none, 
depending on 9. In the latter case, the contact angle is 
too large for forming a liquid surface between the spikes 
and troughs which complies with the boundary condi- 
tion, eq. pj. If, however, A(h) intersects the /i-axis, the 
slopes of the zeros decide upon the stability of the cor- 
responding solutions. This can be seen by appreciating 
that A may be interpreted as a deviation from the force 
balance expressed by eq. Q (35]. For the left zero, which 
is marked by an open circle in the figure, a displacement 
of the three-phase contact line would give rise to an im- 
balance of wetting forces which drives it further away 
from the zero. The opposite is true for the right zero, 
marked by the closed circle. The latter therefore corre- 
sponds to the stable solution, and thus to the adsorbed 
film thickness which will develop. 

All this is in marked contrast to Gaussian roughness, 
for which o\ and 02 are independent of h, with 02/ o\ = 
4/ V V/i [3U]. This is in fact a dramatic difference, as A 
would then just be a horizontal straight line which lies 
either above or below the h axis depending on 9. As a 
consequence, the whole structure we are developing here 
would be absent. 

The graph of A(h) makes a first contact with the hor- 
izontal axis when 9 reaches 



(18) 



At this point, the formerly dry substrate is covered with 
a liquid film of 'thickness' h p — argmax^/cri). A quan- 
tity of particular interest is the total liquid volume, V, 
adsorbed at given 9 and H . This is related to h via 



V= (h- f) P (f)df 



(19) 



and can be evaluated if h and p(f) are known. However, 
we continue here to discuss h instead, since it is more 
accessible through the formalism developed above. 

It is important to note that 9 P always lies above 9w- 
To see this, we note that r = 1/ cos9w ~ 1 + \$w- 
Furthermore, r = (\A + |V/| 2 ). Thus we have 



9 2 



02(h)p(h)dh < max(<72) 



(20) 



On the other hand, (72 > erf , such that 

9 = max — I > max((T2) 



(21) 



From eqs. (20) and (21 ), it follows directly that 9 V > 9\y- 
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FIG. 2: Graphic construction for solving eq. (12 1. The dashed 
line represents the l.h.s. of eq. (1121). 



Let us now consider the system off coexistence, again 
invoking eq. (12) as the condition determining h. A 
graphical solution of eq. (12) is sketched in Fig. [2b. As 



long as 9 > 0, both relevant zeros of A(h)L(h) (represent- 
ing the solutions of eq. (12) for H — 0) lie well within 
the interval [h-,h + ]. For H > 0, the closed circle in- 
dicates again the stable solution. Obviously, the two 
points of intersection will merge when the dashed and 
solid curves touch each other only in a single point. This 
occurs at a certain curvature H p (8) of the liquid sur- 
face. For H > H p , solid and dashed curve meet only for 
h — > h-: there is no liquid adsorbed, and the substrate 
is dry. Hence the average position of the liquid surface, 
h, jumps discontinuously at H — H p . It is clear from the 
construction that H p decreases monotonely with 9. 

As H is reduced below H p , h increases continuously 
until at coexistence it reaches a value corresponding to 
the right zero of AL. Because of the phenomenologi- 
cal similarity of the jump in adsorbed film thickness to 
the prewetting transition encountered in standard wet- 
ting scenarios on flat substrates [5] , it has been proposed 
to term this transition 'Wenzel prewetting' [3(Jj. When 
the microscopic contact angle is varied, a prewetting line 
results, which is shown in Fig. [3] as the solid curve. As in 
the usual prewetting scenario, this line ends in a critical 
end point, when the solid and dashed curves in Fig. [2] 
intersect in only a single point. It is readily appreciated 
from the construction sketched in Fig. [2j however, that 
this can occur only for 9 < 0, which lies outside the phys- 



ically accessible parameter range. 

e 




FIG. 3: The wetting phase diagram wetting on homogeneous, 
isotropically rough surfaces. The most prominent feature is 
the existence of a 'Wenzel prewetting line' (H p (9)), at which 
the adsorbed average him thickness jumps discontinuously 
from zero to a hnite value. 

The variation of h along liquid/vapour coexistence, as 
8 is gradually decreased, can be directly read off Fig. [2^,, 
by inverting the right wing of A(h). This is sketched 
in Fig. [4] as the solid curve which extends between 9w 
and 6 P , and continues in dashed below 9\y- The liquid 
film first appears through a discontinuous jump at 9 P and 
increases gradually as Oyy is approached. As 9 < 9yy, the 
liquid surface configuration which is bound to the surface 
topography through eq. (|3| becomes metastable (dashed 
curve in Fig. RJ, and the global minimum of the total 
free energy corresponds to the 'detached' liquid surface, 
or bulk liquid adsorption (vertical bold line in Fig. p3j). 

It is interesting to discuss the expected shape of ad- 
sorption isotherms, i.e., to consider h(H) at constant 9. 
As it is obvious from Fig. [2b, h will remain finite for all 
positive H, reaching the value indicated in Fig. [4] by the 
sigmoid curve (solid and dashed) for 9 < 9 p and H — >• 0. 
This is remarkable, since there is complete wetting all 
along the coexistence line for 9 < 9w, such that for stan- 
dard wetting scenarios, one would expect the adsorption 
isotherm to diverge continuously as H — > (2]. In the 
present setting, this is only the case if 9 = 0. 

Let us finally get back to the macroscopic contact an- 
gle, m acro- If F(9) is the free energy per unit area of 
the sample for 9 > 9w, force balance at the macroscopic 
contact line yields 



cos 9„ 



F(9) - r 7si 

7 



(22) 



Since F varies continuously with 9, there will be a simple 



zero of 1 — cos 9 r , 



at 6 = 9w- Hence 0„ 
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FIG. 4: Position of the liquid surface (solid curve) and macro- 
scopic contact angle (dash-dotted curve) as a function of the 
microscopic contact angle at liquid-vapour coexistence. As 9 
comes below 9 P , a macroscopic amount of liquid accumulates 
in the troughs and valleys of the roughness. At 9 < 9w, the 
liquid surface detaches completely from the substrate, form- 
ing a bulk liquid layer of arbitrary thickness. The sigmiod 
solid curve corresponds to the right wing of A(ft), displayed 
in Fig. [2k. 



at 9w as \J9 — 9\y, in close analogy to first order wetting 
[2] (dash-dotted curve in Fig. EJ). At 9 = 9 P , where h 
jumps discontinuously to h p , one might at first glance 
expect a jump as well in 9 macro . However, since at 6 P the 
dry substrate is in coexistence with the patchy film with 
'thickness' h = h p , the free energies of the dry surface 
and the wetted surface must be equal, such that F(9) 
is continuous even here. Consequently, there will be no 
jump in 9 macro at 9 p . 

We should not close without mentioning that the pres- 
ence of roughness gives rise to substantial equilibration 
times, as saddle points and extrema occurring at ele- 
vations h Q G [/i_,/i + ] provide effective pinning centers 
[TJ [5J [33J . This will also result in contact angle hystere- 
sis for 9 macro . However, transport through either the gas 
phase or through the molecular adsorbed film [TTJE2] "vvill 
always allow equilibration over manageable times, since 
the lateral distances involved are never larger than the 
lateral length scale of the roughness. 
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